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A series of experiments have been performed to Investigate the oscillatory flow 
pressure and pressure drop through complex geometries. These experiments were 
conducted at the CSU-SLRE facility which is a horizontally opposed, two-piston, 
single-acting engine with a split crankshaft driving mechanism. Flow through a 
rectangular duct, with no insert (obstruction), was studied first Then four different 
inserts were examined: Abrupt, Manifold, Diverging Short and Diverging Long. 
The inserts were mounted in the center of the rectangular duct to represent 
different type of geometries that could be encountered in Stirling machines. The 
pressure and pressure drop of the oscillating flow was studied for: 1) different 
inserts, 2) different phase angle between the two pistons of the engine (zero, 90 lead, 
180 and 90 lag), and 3) for different piston frequencies (5, 10, 15 and 20 Hz). It was 
found that the pressure drop of the oscillatory flow increases with increasing 
Reynolds number. The pressure drop was shown to be mainly due to the gas inertia 
for the case of oscillatory flow through a rectangular duct with no insert. On the 
other hand, for the cases with different inserts into the rectangular duct, the 
pressure drop has three sources: inertia, friction and local losses. The friction 
pressure drop is only a small fraction of the total pressure drop. It was also shown 
that the dimensionless pressure drop decreases with increasing kinetic Reynolds 
number. 
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Cross-sectional area 
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Kinetic Reynolds number 
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Area-averaged fluid velocity 
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Pressure drop factor due to skin 
friction 
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Thickness of the insert 
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Length of the channel 

X/nax 


Swept volume of the piston 
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Displacement of the fluid — — 

4 m 
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Mass of working fluid 
Pressure of fluid 
Pressure drops 

Cycle-averaged pressure drop 

Dimensionless cycle-averaged 
pressure drop 
Radius of the piston 


Greek Letters 
0 

P 

(O 
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Angle shown in Figure 3.24 and 
3.25 

Density of fluid 
Frequency of the piston 
Kinetic viscosity of fluid 


£loc > 4ioc,e > 4loc,c Coefficient of local fluid 
E E E resistance 

J Vtot* 75,c 


L INTRODUCTION 

T HE oscillatory flow pressure drop in the regenerator and heat exchangers plays an important role in 
the Stirling engine performance. The understanding and prediction of the oscillatory flow pressure and 
pressure drop, in engine-like components, are crucial in the optimal design of regenerative thermal 
machines and other reciprocating machines. Currently the modeling of most of the components in these 
regenerative thermal machines relies on correlations derived from steady-state unidirectional flow 
assumptions. Therefore, a more detailed analysis of the heat transfer and pressure drop associated with 
oscillating flow is required for all components of regenerative thermal machines. 

Tanaka, Yamashita and Chisaka 1 described the flow friction and heat transfer characteristics of 
various types of regenerator materials (wire netting, sponge metal and sintered metal) under oscillating 
flow conditions. Zhao and Cheng 2 investigated oscillatory pressure drop characteristics in packed columns 
(composed of three different sizes of woven screen) subjected to a periodically reversing flow of air. Later, 
Lee, Kang and Lee 3 investigated the effects of the regenerator materials on the regenerator effectiveness in 
an oscillating flow. Ju, Jiang and Zhou 4 designed a dynamic experimental apparatus to investigate 
oscillating flow characteristics in a regenerator. 

The Darcy friction factor is related to the pressure drop over a length L by the following equation: 

f APP k 0 ) 

m 

where AP is the pressure drop, u is the cross-sectional mean flow velocity and L is the distance between 
the two positions where the pressure is estimated. 

For steady laminar flow in circular pipes, / is given by: 

/•,« i (2) 

J Re 

For steady turbulent flow regime, / is estimated by: 

, = £3164 (3) 

1 Re #JS 

In a duct with abrupt expansion or contraction in the cross section, local pressure losses occur. The 
local resistance coefficient is defined as: 
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For a uniform velocity distribution over the cross section of the upstream channel, the local resistance 
coefficient of the expansion can be calculated from the Borda-Camot formula (Vennard 12 , 1961) as 
follows: 



For the contraction, the local resistance coefficient (Fried and Idelchik 13 , 1989) is given by: 



For steady flow passing through ducts or pipes with different geometries. Fried and Idelchi 13 published 
“Flow Resistance: A Design Guide for Engineers”, in which practical data for the evaluation of internal 
fluid pressure losses are provided, Su and Lin 5 conducted experimental studies for flows in convergent and 
divergent ducts of rectangular cross section. 

On the other hand, for the oscillatory flow condition, Ni and Gouth 6 ,1997 presented a discussion on 
the dimensionless groups governing oscillatory flow in a baffled tube. Smith and Swift 7 (2003) performed 
some experiments on oscillatory flow at the abrupt transition between a two-dimensional channel and 
essentially infinite space. 

Chen and Griffin 8 (1983) studied the pressure drop correlations for Stirling engine applications. Their 
study indicated that the regenerator pressure drop has a significant effect on the engine efficiency and 
power output. And the pressure drop in the regenerator is nearly a linear function of velocity. Also, they 
found the ratio of friction factor between oscillatory and steady flow only varies with the kinetic Reynolds 
number. 

In this paper we describe some of the research effort that was conducted at Cleveland State University 
(CSU) with the following objectives: (1) to design and build a data acquisition system for the CSU-SLRE 
(Stirling Laboratory Research Engine) with the capability to measure the pressure, pressure drop and crank 
angle with high speed (up to 60 Hz), (2) to analyze the pressure and pressure drop signals in a rectangular 
channel for both no-insert and different inserts (of different geometries), for various Reynolds numbers and 
for different phase angles between the two pistons (see Figure 1), (3) to study the variation of 
dimensionless pressure drop as a function of kinetic Reynolds number for different inserts and (4) to 
compare the local flow resistance loss i.e. “minor” or “end effect” loss, of oscillating flow with steady flow 
at the same geometry. 


IL EXPERIMENTAL DETAILS 

The experimental facility is shown in Figure 1 . The SLRE, as originally designed and fabricated, is a 
horizontally opposed, two-piston, single-acting engine with a split crankshaft driving mechanism. The test 
section is a parallelepiped with 17.3 mm width, 173 mm height and 300 mm length (see Figure 2.), In these 
tests, the working gas (air) passes through the test section mounted between two transition sections, which 
in turn are connected to the two cylinders (with two pistons) of the engine. The transition sections, 
connecting the piston head to this test section, have a rectangular cross section (17.3 mm width and 173 
mm height) on one end and a square section (52.9x52.9 mm) on the other. This design was made to ensure 
smooth transition of the flow from piston head to the test section. Thus the maximum volume of the test 
section is 0.002687 ra 3 with both of the two pistons at BDC, while the minimum volume is 0,002237 m 3 
with both of the two pistons at TDC. 
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Pressure drop transducer 



Figure 1. Schematic Diagram of the Experimental Apparatus. 



Figure 2. Rectangular Duct (no insert). 



Figure 3. Rectangular Duct with an “Abrupt” insert. 



Figure 4. Rectangular Duct with a “Manifold” 


4 

American Institute of Aeronautics and Astronautics 





I 

1.730 - 




^//////////////////^ 


L 7.440 1 

25.400 — 


Figure 5. Rectangular Duct with a “Diverging” short 
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Figure 6. Rectangular Duct with a “Diverging” long insert. 


Figure 2 shows the rectangular test section while the different type inserts used in this study are shown 
in Figures from 3 to 6. These inserts were mounted in the central part of the rectangular test section. 
Oscillatory flow pressure and pressure drop through the rectangular duct (with no insert) were studied first. 
Then four different inserts were used: Abrupt (ABRU) in Figure 3, Manifold (MANI) in Figure 4, 
Diverging Short (DIVS) in Figure 5 and Diverging Long (DIVL) in Figure 6. Note that insert dimensions 
are defined in these figures. The pressure and pressure drop of the oscillating flow was studied for: 1) 
different inserts, 2) different phase angles between the two pistons of the engine (zero, 90 lead, 180 and 90 
lag), and 3) for different piston frequencies (5, 10, 15 and 20 Hz). The zero degree phase angle difference 
between the two pistons means that one piston is at the Top Dead Center (TDC) when the other piston is 
also at the TDC. The 180 degree means that one piston is at the TDC when the other is at the Bottom Dead 
Center (BDC). After one piston was chosen as the reference piston, the 90 degree lead or lag, represented 
the phase of the 2 nd piston’s motion relative to that of the reference piston. 

Different experimental cases were represented with a particular code such as: insert-frequency-phase 
angle difference, e.g. [ABRU-5-90 lead] indicates an abrupt insert, frequency = 5 Hz and the phase angle 
difference = 90 lead. 


IIL RESULTS AND DISCUSSION 

In this section, we present the experimental results of pressure and pressure drop through the 
rectangular duct with no insert and then with different insert (obstruction) geometries. 

A, Pressure data for 180 degree phase angle difference with no inserts 

The 180 degree phase difference cases (when one piston is at TDC, the other will be at BDC) have the 
lowest pressure change inside the test section, since for this case the two pistons are moving in the same 
direction with the same speed. The measured pressure versus crank angle data for the no-insert (N-5, 10, 
15, 20-180) cases are presented in Figure 7. It is shown that the pressure data for different frequencies with 
no-insert are very similar. Figure 7 shows two peaks in one cycle (rather than one) and they do not have the 
same magnitude. These two peaks are due to the fact that the piston motions are not perfect sinusoidal 
functions. When the two pistons move at “180 phase angle” difference, the total volume of the working gas 
changes slightly. Also as seen from this figure, the signals are not symmetric about the 1 80 crank angle. 
This is due to the fact that the phase angle difference between the two pistons is not exactly 180 degrees. A 
small error in the phase angle, such as phase drift from the belts, will cause such an unsymmetrical result. 
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Figure 7. Experimental Pressure Signal of N-5, 10, 15, 20-180 Cases 


Using the motions of the two pistons to calculate the pressure inside the test section, Figure 8 shows 
the results of the pressure variation inside the test section for 1 80 degree phase angle difference between 
the two pistons. On the other hand. Figure 9 shows the pressure variation with 182 degree phase angle 
difference between the two pistons. The two degrees phase angle difference resulted in a change of the 
negative pressure “peak” from being 600 Pa to two peaks at 400 and 800 Pa. Further more, there are 
differences among the four pressure signals (frequency = 5, 10, 15 and 20 Hz) as shown in Figure 7. This is 
because the pressure transducer was not installed exactly in the middle of the test section. As the frequency 
(flow velocity) goes higher, the pressure difference between the middle of the test section and the location 
of the pressure sensor will increase. The biggest difference appears at the TDC and BDC, because the 
maximum inertial pressure drop happens at TDC and BDC. 



Figure 8. Calculation of Pressure Signal at Figure 9. Calculation of Pressure Signal at 

1 80 Degree Phase Angle Difference 1 82 Degree Phase Angle Difference 
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B. Pressure drop of oscillatory flow through rectangular channel (test section without insert) 

Figure 10 shows the pressure drop versus crank angle for the cases of no insert, 180 phase angle 
difference and frequencies 5, 10 15 and 20. For our case we have a channel with Dh (hydraulic diameter) = 
0.0309 m, L (length in the flow direction) = 0.3 m and aspect ratio = 10, or an L/Dh value of 9.7. 
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Figure 10. Experimental Pressure Drop Signal of N-5, 10, 15, 20-180 Cases. 


For steady flow we used two correlations for friction factor: (1) fully developed flow and (2) 
developing flow. The correlation used for the steady flow condition, fully developed and with aspect ratio 
=10 is: 

/•Re = 21.17 (7) 

For the developing laminar flow condition, the correlations between the L/(D h *Re) (from 0.002 to 
0.1) and f*Re were shown schematically for different aspect ratios in Figure 41 in the book “Laminar Flow 
Forced Convection in Ducts”, by Shah and London 14 . From their results: 

/ • Re a 65 (8) 

Table 1 shows the maximum pressure drop obtained experimentally together with friction pressure 
drop obtained from steady state correlations. In comparing between the oscillatory and steady flow cases 
we used the quantities obtained at the maximum Reynolds number. Taylor and Aghili 9 (1984) studied the 
friction factor correlations between the steady flow and oscillatory flow in tubes. They found the oscillatory 
flow friction factor was 3 to 4 times larger than the steady flow when L/D equals to 13.7 and Reynolds 
number is from 10 3 to 10 5 . We can see from Table 1 that the maximum experimental pressure drop, at the 
maximum oscillatory-flow Reynolds number were hundreds of times larger than the unidirectional-flow- 
friction pressure drops. Therefore, we conclude that the oscillatory-flow pressure drop measured in the 
rectangular duct used in this study was not a frictional pressure drop. 

The velocity and acceleration profiles of each piston can be calculated from its displacement. By 
assuming that the working gas moves at the same pace as the pistons (for the real situation, there is some 
lag between the motions of the gas and the pistons), the gas inertial pressure drop can be calculated from 
the known acceleration by: 

AP = p'l' a (9) 
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Table 1 The maximum experimental pressure drop and the calculated steady-flow 
frictional pressure drop, both corresponding to the maximum oscillatory flow 
Reynolds number, at different frequencies. 


CO 

(Hz) 

Re, max 

Fully developed flow 
AP (pa) - Steady Flow 

Developing flow AP 
(pa) - Steady Flow 

Measured max AP 
(present work) (pa) 

5 

3319 

0.088 

0.269 

16.95 

10 

6637 

0.175 


60.49 

15 

9956 

0.263 


137.14 

20 

13275 

0.350 


239.28 



Figure 1 1. N- 15- 180 Pressure Drop Comparisons between No-Insert (N- 15- 180) Experimental Data and the 

Inertial Pressnre.-Dron Calculation Results 

Figure 11 compares the pressure drop curve (N-15-180) obtained in these experiments with the one 
calculated for the inertial pressure drop using equation (9). As can be seen, although there’s a difference in 
the magnitude (maximum difference 24%) and a slight difference in the phase, these two curves are similar. 
The difference in the magnitude is thought to be due to the fact that a fraction of the measured pressure 
drop is due to skin friction; the mismatch in the phase is, then likely, due to the fact that the maximum 
value of the frictional pressure drop occurs at 90 and 270 degrees phase angle (at maximum piston velocity) 
while the maximum value of the inertial pressure drop occurs at 0 and 180 degrees phase angle (at the 
maximum piston acceleration). The overall good agreement shown in Figure 11 further justifies the 
conclusion that the pressure drop of oscillating flow in the rectangular duct under investigation, is mostly 
due to the gas inertia. 

C. Pressure drop of oscillatory flow through rectangular channel with different inserts 

Figures 12 to 15 show the experimental oscillatory pressure drop for the four different inserts 
inside the rectangular duct. The pressure drops for oscillatory flow through the rectangular duct with no 
insert were mainly from the inertia as shown before. However the pressure drop for oscillatory flow 
through the rectangular duct with inserts will be shown to come from three different sources: inertia, local 
losses (minor or end effects loss) and skin friction (which is thought to be small based on the previous 
comparison for no insert). 
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Figure 12. Experimental Pressure Drop Results of Figure 13. Experimental Pressure Drop Results of 
Four Abrupt Inserts (ABRU-5, 10, 15, 20- Four Manifold Inserts (MANI-5, 10, 15, 

1 80) Cases. 20- 1 80) Cases. 
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Figure 14. Experimental Pressure Drop Results 
of Four Short Diverging Inserts (DIVS- 
5, 10, 15, 20-180) Cases. 
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Figure 15. Experimental Pressure Drop Results 
of Four Long Diverging Inserts (DIVL- 
5, 10, 15, 20-180) Cases. 


For the oscillatory flow condition, the pressure drop from the skin friction is just a small fraction of the 
total pressure drop in the rectangular duct. For inertial pressure drop, the maximum signals appear at the 
phase angle of 0 and 180 (at the maximum piston acceleration). For the pressure drop due to local loss, the 
maximum signals appear at the phase angle of 90 and 270 (at the maximum piston velocity). Since the 
actual situation is that the pressure drop is due to the total effect of the inertia and the local loss in this 
series of experiments, the phase angle of the maximum magnitude of the signal is located between 0 and 90 
and between 180 and 270, as shown from Figures 12 to 15. 

A comparison of the pressure drop among different phase angle differences between the two pistons 
(90 lead, 180 and 90 lag) is shown in Figure 16 for the abrupt insert and frequency = 20 Hz. At the same 
frequency, the maximum pressure drop over the cycle for the 90 degree lead or lag, is shown to be less than 
the maximum pressure drop for the case of the two pistons moving at 180 degree phase difference (for 
which the working gas experience the maximum velocity and acceleration of all phase differences). 
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Figure 16. Pressure Drops of Different Phase Difference between the Two Pistons 
with Abrupt Insert. 


D. Dimensionless Cycle-Averaged Pressure Drop 

To get a relationship between the dimensionless cycle-averaged pressure drops and the kinetic 
Reynolds number, the cycle-averaged pressure drop was calculated using the following equation (It should 
be noted that this pressure drop was taken at 1024 equally space crank angle position over the cycle): 

1024 

AP = £|A/>| (10) 

n — 1 

The dimensionless mean pressure drop factor and the kinetic Reynolds number are defined as follows: 


AP' = ^ 

PLUU 2 

(11) 

IncoDl 
Re„= 

(12) 


where co is the frequency of the piston: 5, 10, 15 and 20 Hz. L is the distance between the two points of 
pressure drop measurement. 

And, 


4 wb 

h 2{w + b) 

(13) 

rr 2 ^max 

max 2 

(14) 

'^max A 

A m 

tn 

(15) 

Vpis S 

(16) 

A =(w-h)-b 

(17) 


where r is the radius of the piston, s is the stroke of the piston. W and b are the width and the height of the 
rectangular duct respectively. The thickness of the insert is denoted as h. The relationship between the 
dimensionless cycle-averaged pressure drops, with inserts and the kinetic Reynolds numbers are shown in 
Figure 17. This figure shows the dimensionless pressure drop decreases with increasing kinetic Reynolds 
number for the cases of oscillatory flow passing through rectangular duct with inserts. At first this may 
seem a bit surprising. However a look at the definition of dimensionless pressure drop in equation (11), 
indicates that the pressure drop is divided by the square of piston velocity and piston velocity is 
proportional to the frequency. 
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Figure 17. Relationships Between the Cycle- Averaged Dimensionless Pressure Drop and Kinetic 
Reynolds Number. 


E. Comparison of local pressure drop between oscillatory and steady flows 

In order to compare the oscillatory-flow local pressure drop with local pressure drop for the steady 
flow condition, the inertial pressure drop needs to be estimated so that the local pressure drops can be 
segregated from the total pressure drop. 



— ► M,-d 



M^-a M 3 -a 



Figure 18. Relationships Between the Pressure Drop and the Acceleration with Insert. 


For the test section with the abrupt insert, as shown in Figure 18, the pressure drop from the inertia can 
be calculated by: 

(P x - P 2 ) A = M x a + M 3 a + M, a (18) 



M 3 -a 


M, • a M 2 - a 

Figure 19. Relationships Between the Pressure Drop and the Acceleration without Insert 


For the test section with no insert, as shown in Figure 19, the pressure drop from the inertia is: 

(P l -P 2 )A = M l a + M 3 a+M 2 a (19) 
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Since the mean pressure for the above two cases are equal, if M 2 
time of a, that is: f 

M 2 a ' = M 2 a 


is 1/n time of M 2 , a' will be n 


( 20 ) 


The same also applies to the other kinds of inserts. Therefore, inertial pressure depends only on the gas 
acceleration (and therefore the frequency of the pistons) and the mass of the gas being accelerated. 

Figures 12 to 15 show the experimental pressure drop signals for different inserts at different 
frequencies, with the two pistons moving at 180 phase angle difference. Since the pressure drop from the 
inertia only depends on the frequency, the inertial pressure drop for the case of the abrupt insert is 
approximately the same as the pressure drop for the no insert case. Thus, by subtracting the no insert 
pressure drop from the total pressure drop of the abrupt insert case (at the same frequency), the pressure 
drop left is the local pressure drop due to the expansion and contraction inside the test section, as shown in 
Figures 20 to 23. 



Figure 20. Oscillatory Flow Local Pressure 
Losses of Abrupt Insert (ABRU-5, 10, 15, 
20-180 Cases). 


Figure 21. Oscillatory Flow Local Pressure 
Losses of Manifold Insert (MANI-5, 10, 
15, 20-180 Cases). 



Figure 22. Oscillatory Flow Local Pressure 
Losses of the Short Diverging Insert 
(DIVS-5, 10, 15, 20-180 Cases). 


Figure 23. Oscillatory Flow Local Pressure 
Losses of the Long Diverging Insert 
(DIVL-5, 10, 15, 20-180 Cases). 
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Under the steady flow condition, with the abrupt insert, the local resistance coefficient can be 
calculated for different geometries of expansion and contraction using equations 5 and 6, where Di, D 2 and 
u are shown in Figure 24. Di, D 2 are hydraulic diameters of the small section and large section respectively 
and u is the flow velocity in the small section. 


T 


d 2 


< ► U 


Di 


Figure 24. Schematic Diagram of Flow with an Abrupt Change of the Cross-sectional Area. 


To calculate the steady-flow pressure drop in the rectangular duct with manifold insert, the total local 
resistance coefficient was calculated by: 

£,,=2(# fec ,+4j (21) 

For the diverging insert cases, the steady flow local resistance coefficient should be calculated based 
on the smooth transition of the section area, as shown in Figure 25. An analytical formula to calculate the 
steady-flow local resistance pressure drop is only available for the geometry shown in Figure 26. These two 
conditions are similar because the ratios of the contraction or expansion of the cross-sectional area are the 
same for the same 8 in these two figures. Therefore, it is reasonable to use the geometry in Figure 26 to 
approximate the geometry in Figure 25. 



Figure 25. Schematic Diagram of Figure 26. Schematic Diagram of Flow with a 

Flow with a Smooth Change of the Symmetric Smooth Change of the Cross- 

Cross-section Area. section Area. 


'i-A v 


For the geometry shown in the Figure 26, under expansion condition: 
where 

£ =0.01 10' 22 , 7.5° <9 < 35° 

and 0 is in degrees. 

Under the contraction condition, the resistance coefficient is: 


( 22 ) 


l, = 


-0.0125^ j +0.0224^ j -0.00723^- j +0.00444 ^--0.00745 x(f/ r 3 -Iny) -\0y/ r ) 


(23) 


where 

y/ r = 0.017450 (8 in degrees) 

Therefore, for the steady flow passing through diverging insert, the total local resistance coefficient is: 

C 2 ^24) 

where 
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( 25 ) 

4,,2=Ac+4c* (26) 

Figure 27 shows the calculated pressure drop for the steady flow and the cycle-averaged local pressure 
losses from the experiments for different inserts. It should be noted that we used the cycle-averaged 
Reynolds number in this case. Typically the maximum Reynolds would be used (in channel flows and as 
shown earlier in the no-insert cases). However due to the higher harmonics present in the complex 
geometries (i.e. with inserts, see Figure 16) we decided to use the cycle-averaged Reynolds number. The 
local (minor) pressure drops of steady flow are almost the same as the cycle-averaged oscillatory pressure 
drop for the same flow velocity in the rectangular duct with an abrupt insert. The manifold insert cases 
show the greatest difference between the calculation and the experimental result. It is believed that is 
because the distance between the two parts of the manifold insert is not far enough to avoid interference 
between the two parts. The magnitude of the local pressure drops for the diverging short and diverging long 
inserts are very similar. This indicates that the ratio of the expansion (contraction) is the most important 
issue in determining the local pressure drops for the smooth transitions of cross-sectional area. In both 
cases, the experimental results are 1 .2 ~ 1 .3 times the calculation results for the steady flow. 
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Figure 27. Comparison of Local Pressure Drop Between Experiments (Oscillatory Flow) and Calculation 
(Steady Flow) for Different Inserts, for the same Cycle- Averaged Reynolds Number. 



IV. CONCLUSIONS 

The maximum pressure drop and the cycle-averaged pressure drop of oscillatory flow increase with 
increasing piston frequency (and accordingly maximum Reynolds number or kinetic Reynolds number). 
The pressure drop for the oscillating flow in a no-insert rectangular duct is mainly caused by inertia. The 
oscillating-flow pressure drop in a duct with an insert is due to the total effect of the local (or end effect) 
pressure drop and the inertia. For the same frequency and geometry, the case of the two pistons moving at 
the same phase (180 degree) has the highest cycle-averaged pressure drop. The dimensionless pressure 
drops decrease with the increasing kinetic Reynolds number, which seems reasonable when it is noted that 
the denominator of the dimensionless pressure drop is proportional to the square of the frequency. For 
oscillatory flow through a rectangular duct with the inserts considered in this study (abrupt, manifold, 
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diverging), the cycle-averaged local pressure drop under the oscillatory-flow condition is almost the same 
as that under the steady flow condition, for the same geometry. 
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